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Temperature 


Modeling of Directional Solidification 


Diffusion controlled growth 

Frozen temperature approximation T(z)=T 0 +G(z-V 0 t) 

Local equilibrium at interface 

One dimensional case - planar interface growth 



G - thermal gradient in liquid phase 
V 0 - furnace translation rate 
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Governing Equations 


Solute diffusion in liquid phase in 
the coordinate system co-moving 
with the interface 

Interface boundary condition 


Interface velocity 
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Initial Conditions 


At time t=0: 

Interface velocity is zero 

Solute concentration in liquid phase is C 0 everywhere 
Solute concentration in solid is kC 0 
Temperature at the interface position z 0 is T 0 -Gz 0 

Initial transient during directional solidification was treated by Tiller in 1953. Some 
recent works include: 

• W. Huang et al, J.Cryst. Growth 182 (1997) 212-218. 

•D. Ma et al, J. Cryst.Growth 169 (1996) 170-174. 

®A. Karma, A. Sarkissian, Phys.Rev. E 47,(1993) 513-533. 

• B. Caroli et al, J.Cryst. Growth 132 (1993) 377-388. 

• Majchrzak et al, J. Mater. Proc. Techn. 78 (1998) 122-127. 

•Ch. Charach et al, Phys. Rev. E 54 (1996) 588-598. 

« M. Conti, Phys.Rev. E 60 (1999) 1913-1920. 
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Interface Velocity 


Numerics 


Numerical solution by COMSOL 3.3 and FLEXPDE 5 




M. Conti, Phys Rev. E (1999) 
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Volterra Integral Equation 
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Here c 0 (t) it the interface concentration, z(t) is the interface position , 

V(t) is the interface velocity, CJs the concentration in the liquid far from interface 
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Analytic Analysis 


Tiller's case of a constant 
growth rate V=V 0 
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Laplace transform technique 
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Inverse Laplace recovers 
Tiller's formula 



]-b-b erfc ( b\ft j 



Small time asymptotics 
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Interface velocity [cm/sec] 


Results 



Time [sec] 





The 21st Conference on Crystal Growth and Epitaxy-West 
June 8 - 11, 2008, Stanford Sierra Camp, Fallen Leaf Lake, CA 


UAH 9 

M Iftfrarafty cl Alabama In Huntsville 



Conclusions 


• Strong spikes in growth velocity are obtained for the simplest, 
physically sound, start-up model of directional solidification. 

• Spikes are sensitive to the flow in the system and to the latent heat. 

• Integral equation approach can be managed numerically easier than 
the Finite Elements/Volume methods. 

• The spike phenomenon is not uniquely related to the phase field 
model of the solidification. 
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